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Abstract 
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I. INTRODUCTION 



In the calculation of some molecular property, the magnetic shielding for example, one 
usually assumes fixed position of nuclei, the so called Born-Oppenheimer (BO) approxima- 
tion and at the final stage, averages over the appropriate vibration-rotational function. The 
principal question we address in this work is, what are the finite nuclear mass corrections 
to various physical properties such as relativistic energies, polarizabilities or the mentioned 
magnetic shielding evaluated in the BO approximation. The direct nonadiabatic calculations 
are possible only for small molecules and only for simple properties such as the electric dipole 
polarizability l|, |2j . This approach however is not universal and has not yet been applied 
to more complex molecules such as H 2 0, or to the evaluation of the nuclear spin- rotation 
and shielding constants. In this work we demonstrate the applicability of the nonadiabatic 
perturbation theory [3| to obtain in a systematic way the formulae for various physical 
properties of a diatomic molecule, with possible extensions to larger molecules. We rederive 
the known result for the electric static polarizability, the rotational magnetic moment, the 
spin-rotation constant, and obtain the leading finite nuclear mass corrections (which we will 
call in this work nonadiabatic corrections) to relativistic rovibrational energies, the tran- 
sition electric dipole moment, the electric static polarizability, and the magnetic shielding 
constant. These finite nuclear mass corrections are important for comparison between ac- 
curate measurement and precise calculations, for example in the dissociation energy of H 2 
0, 0], the transition electric-dipole moment of LiH jf|, or in the shielding constant of H2 
and isotopomers Q]. 

We demonstrate in this work that the leading finite nuclear mass corrections can be 
conveniently calculated for a fixed position of nuclei, as in the BO approximation, and 
averaged out over the rovibrational function. We do not consider here the second order 
corrections in the nuclear Hamiltonian H Q , but for consistency regard them also as the 
nonadiabatic corrections, although of higher order. 

In Sec. II we define the reference frame and split the nonrelativistic Hamiltonian into 
the electronic and nuclear parts. In Sec. Ill we briefly present nonadiabatic perturbation 
theory on the basis of Ref. include nonadiabatic corrections to the BO wave function 
and derive general formulae for the first and the second order matrix elements. In Sec. 
IV, as a first example, we derive formulae for relativistic recoil corrections to rovibrational 
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energies in diatomic molecules. Other examples, the finite nuclear mass corrections to electric 
properties of molecules are derived in Sec. V, and to magnetic properties in Sec. VI. Finally, 
the obtained results are briefly summarized in Sec. VII. 



II. NONRELATIVISTIC HAMILTON! AN 

We consider a neutral diatomic molecule with the Hamiltonian 

„2m 2rriA 2m B 

where the summation index a goes over all electrons and A and B refers to nuclei. In order 
to derive formulae for nonadiabatic effects, one must fix the reference frame. We start with 
the laboratory frame {Ra, Rb, r a }, and change variables to {R, Rg, x a } according to 

Ra = Rg + £b R, (2) 
Rb = Rg — £a R, (3) 
f a = R G + x a , (4) 

with the relative position of nuclei R = Ra~Rb, and the arbitrarily chosen on the molecular 
symmetry axis new frame origin R G = c^-Ra + £bRb, where €a + £b = 1- The conjugate 
momenta are related by 

Pa = € A P G + P-€ A ^q a , (5) 

a 

Pb = e B P G - P - e B Qa, (6) 

a 

Pa = q a , (7) 

— * — » — » 

where P = —i Vr and q a = —i V ' Xa . The nonrelativistic wave function with the vanishing 

— * — * — * 

total momentum does not depend on R G , thus <p = <f>(x a , R) and since P G commutes with 
H when expressed in new variables, it can be set to 0. Hamiltonian H in new variables 
becomes 

H ~ ^2~^ + V + {2mA + 2m B ) + {2mA + 2mJ ^ a 



a 



( £a _£b_\ 

\rriA m B J 



(8) 
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The last term in the above is transformed by the unitary transformation 

H = e - 4v He i ' p = H-i[(p, H] + ... (9) 

where 

^ = ™(— -— )E*VP, (io) 

\m A m B J V 

and 

e~ iv £e i * = £ + m(^--^) P\ (11) 
\m A m B J 

with higher order 0(m/M) 2 terms in the electron nuclear mass ratio being neglected and 
M is equal to m A or m B . As a result of this transformation Hamiltonian takes the form 

H = H cl + H n , (12) 

ff - = ?^ + V « (13) 

-m(^-^-\ J2x a -V R (V) (14) 
\m A m B J a 

= K + K, (is) 

where if£ includes the first term and H£ the two remaining terms. This form of the nuclear 
Hamiltonian is convenient for the calculation of nonadiabatic effects. Moreover the freedom 
in choosing e A ,B will be used in order to simplify formulae for nonadiabatic corrections to 
electric and magnetic properties. Later, we will need the angular momentum operator J, 
which for states with the vanishing total momentum is defined by 

J = ^2(f a - Rcm) XPa+ {Ra ~ Rcm) XPA + {R-B ~ Rcm) X p B , (16) 
a 

where Rcm is the molecular mass center. In new variables the operator J 

J = ^2x a xq a + RxP = J cl + J n (17) 

a 

is split into electronic J c \ and nuclear J n parts, and this J is not being modified by the 
unitary transformation of Eq. (flOl) . 
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III. NONADIABATIC PERTURBATION THEORY 



The total nonrelativistic wave function <p of an arbitrary molecule is the solution of the 
stationary Schrodinger equation 

[H-E]\<j>)=Q, (18) 

with the Hamiltonian H being a sum of the electronic H e \ and nuclear H n parts, Eq. (1121) . 
In the adiabatic approximation = a , where 

<P^x,R) = <p cX {x) X {R) (19) 

is represented as a product of the electronic wave function <p e \ and the nuclear wave function 
X- We note, that <j) e \ depends implicitly on the nuclear relative coordinate R. The electronic 
wave function obeys the clamped nuclei electronic Schrodinger equation 

[# el -£ e i(#)]|0 e i)=O, (20) 

while the nuclear wave function is a solution to the Schrodinger equation in the effective 
potential generated by electrons 

[H a + S a (R) + S cl (R) - E a ] \ X ) = 0, (21) 

where 

£ a OR) = (0el|#n|0el) el . (22) 

In the nonadiabatic perturbation theory, the total wave function 

= a + <50 na = 0el X + <50na (23) 



is the sum of the adiabatic solution and a nonadiabatic correction. The nonadiabatic cor- 
rection <50 na is decomposed into two parts 

<$0na = 4>c\ SX + <5'0na , (24) 

obeying the following orthogonality conditions 

(5>na|0el)el = 0, (25) 

(5x\x) = 0, (26) 
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which imply the normalization condition ((/> a \<f>) = 1- 

In the first order in H n of the nonadiabatic perturbation theory one has 

1 



|<5'0na) 



(1) 



and in the second order 

|5V>na> (2) 



1 



(£ e l — H e \) 



[H n + £ el — E a 



-H n \cj) e} x), 



-H D |0 el x), 



(27) 



(28) 



(S e i — H e \)' (£ C \ — H C \) 

where l/(£ e i — H e {)' denotes resolvent with the reference state e i subtracted out. The total 
nuclear function x + $X satisfies the effective Schrodinger equation which includes adiabatic 
and nonadiabatic corrections 3j. Thus the nonadiabatic wave function can be recovered 
order by order in the perturbative approach. Eqs. fl27|28l) involve H n , and thus derivatives 
with respect to R. These derivatives can be calculated with the help of the following formulae 

1 



V 



R 



Vfl|&a) 
1 



(£ el — H { 



el J 



(£ i — H c \) 

1 



1 



7 V H (V)|0 el ), 



(29) 



1 



(£ c l — -"el 
-|0el)(0el|V fl (l/) 



(£ C 1 — H c \) 

v«(y)|0 el )(0 e i| 
i 



(30) 



It has been shown recently in [3|, that the application of these formulae allows for a significant 
improvement in the numerical accuracy of adiabatic and nonadiabatic corrections in H2 
molecule, and this probably will hold for any diatomic molecule. Alternatively, one may use 
the formula 

V R = n(n- V fl ) -n x (n x V R ), (31) 

where n = R/R, with only the first radial part in above equation replaced in terms of the 
derivative d(V — £ e \)/dR, namely 



'R\<Pe\ 



dV , 



n 



(£ cl -H cl )' dW^ 1 B 
For example the adiabatic correction to energy becomes 



) - —nx J n |0 e l). 



(32) 



^>el|#n|0el) 
UK ' 



J 2 



2m n R 2 



\4>c\) + (<Pe\ 



dV 



1 



dV, 



dR (S cl -H el )' 2 dR> 



1). 



(33) 



where m n is the nuclear reduced mass and for £ electronic state J n |0 e i) can be replaced by 
— J c \\4> e \). In this way one avoids summation over intermediate states with the II symmetry. 



A. First order matrix elements 



We will use here this nonadiabatic perturbation theory to derive the finite nuclear mass 
corrections to various matrix elements in the general form. Later we will analyze specific 
examples. Consider the Hermitian electronic operator Q (no derivatives with respect to 
nuclear variables, for example the relativistic correction to kinetic energy of electrons) and 
its matrix element between (different) rovibrational states. In the BO approximation this 
matrix element can be represented in terms of the electronic matrix element nested in the 
nuclear matrix element, namely 



= (0eiXf|Q|0ciXi) = (xfl(Q)^IXi), 

(Q)<® = (Q) el = (0 cl |Q|0 el >. 



(34) 
(35) 



We will show that the same holds for nonadiabatic corrections to this matrix element, which 
are 



(Q) (1) = <0clXf|#n- 



1 



77 <5 1 0ci Xi) + (<j>dXf\Q 



1 



-H^xi). (36) 



Although in this work we consider only the first order 0(m/M) corrections, let us present 
here the second order corrections to the diagonal matrix element to demonstrate the appli- 
cation of the nonadiabatic perturbation theory 



(Q) 



(2) 



Xl#n 



1 



7j (H n + £ c \ — E a 



yQ\<Pe\ X) 



(£ ei — H c i)' (£ c \ — H c \j 

+<^i X\Q ,~ 1 „ v ^ + S * ~ E ^ (P 1 n v H ^ 
(£ el — n e i) (£ el — n e i) 

+(<f>ei x\H n - c (Q - (Q> (0) ) H n \<f> el x). 



(37) 



(£ el — H c \)' (£ c i — H e ij 

Additional corrections due to 5\ i n Eq.( l24l) can easily be included in the (Q)^ and (Q)^ 
by replacing x by X + $X an d will not be considered any further. Let us return now to the 
leading order correction to the matrix element in Eq. (1361) 

1 „ „ 1 



(Q) 



(i) 



J d?R{{xlXi 
V(Xf Xi 



(£ el — H c 



-Q\cf> el ) + (cj) el \Q 



(£ el — H c 



w \H n (pci) 



2 m n 

(xiVx? -Xi Vxi 
2 m„ 



\t c \ — ti c \)' {t e i — li e i) 



(38) 



(V 



R<Pel 



(£ el — /J, 



el J 



-QM-(MQ 



(£cl - #el)' 



IV 



R<Pel 
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and consider two special cases. If Q is a real operator, then the third term vanishes and 
with the help of integration by parts we obtain 



(Q) {1) = (MQ)i ] \xi), 



(Q) 



a) 

el 



(£ c i — H e \) 



7 Q|0el> + (MQ 



1 



(£el - #el)' 



H n |0el), 



(39) 
(40) 



where for arbitrary vfj e \ and ip' el : 



= (V«^i|VH^ci)/(2m n ) + (^il^l^ei)- 



(41) 



This case of the real Hermitian Q finds applications in studying relativistic corrections to 
rovibrational energies and to all electric properties. If Q = Q is an imaginary vector operator 
composed of electronic operators and R, such that R- Q = 0, then the first two terms in Eq. 
(158]) vanish and 

1 r d 3 R 



(Q 



n(i) 



2m r 

x 



R 2 



[G£xf)*Xi + x? (■£*)] 



>el|Q J 



(£ c i — H e \) 



7 | j^0 e i) + (JIM 



(£el — -f^c 



(42) 



where J n is defined in Eq. (fl7|) . Let us assume that e i is a £ state, then J n |0 el ) = — J e \ \4> e \) 
and 



Xf 



J' 



1 



(43) 



m n i? 2 ^ Cl (£el-^el 

This case of imaginary Hermitian Q finds application in studying the magnetic properties. 



B. Second order matrix elements 



Consider the second order matrix element with two arbitrary electronic operators Q\ and 
Q 2 , Let us assume that (0 e i|Qi|0ei) — and introduce the notation 



(Qi Q 2 > 



l 



In the leading order of the nonadiabatic perturbation theory this matrix element is 



(44) 



(Qx Q 2 > 



(0) 



>el\Ql 



(£ el — H e 



w^l^el) +C.C, 



(45) 

(46) 
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and the nonadiabatic correction is 

(QiQ 2 > (1) = (xMQi (f . 1 „ ,, {H n + 8 cX -E & ) L— Q 2 |0 dX ) 
+ (X 0oi|^n 77; — ^-JTy Qijc — ~7Tv X) 

{tel — rlel) [^el — Mel) 

HxMQi (P 1 f T,, Q2 (P —-y, H a \(f> el x) + c.c. . (47) 

This correction can also be rewritten in terms of the nested electronic matrix element, 
namely 

(QiQ 2 ) (1) = (xliQiQ^lx) (48) 
(Qi Q 2 )S ) = (0eilQi ^ 1 (H n -ej 1 q 2 |&i) 

(Al — -"el J l^el — -"el J 

+ (0el| An 77 —fTyQlJ^ ^rT7<52|0el) 

[C e i — ilelj (£-el ~ -"clj 

+ 77J -TTv * , H n |0 cl ) + c.c. . (49) 



These formulae will be used in the calculations of the nonadiabatic corrections to the shield- 
ing constant. The more general case with {4>ei\Qi\4>ei) 7^ for Hermitian real operators Qi is 
considered in the following section using a slightly different approach. 

C. Diagonal matrix elements with real Hermitian operators 

The finite nuclear mass corrections to the diagonal matrix element of a Hermitian and 
real operator Q can obtained by taking a derivative Sq with respect to Q of the nuclear 
Schrodinger equation, which includes the diagonal adiabatic correction S a in Eq. (12ip 



6 Q [H n + £ a (R) + S cl (R) - E a ] |x> = , (50) 

that is 

[H a + £ & {R) + S cl (R) - E a ] \5 qX ) + [S Q S a (R) + 5 Q S cl (R) - 5 Q E a ] \ X ) = . (51) 

Taking the product with (x\ on the left hand side, one obtains the matrix element with the 
leading finite nuclear mass corrections 

(0|Q|0> = S Q E « 5 Q E & = ( X \6 Q £ e i(R) + S Q £ a (R)\ X ). (52) 
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The perturbation of electronic energies £ e \(R) and £&(R) due to some operator Q can be 
obtained using the standard Rayleigh-Schrodinger perturbation theory and the result 

(m<P) = (x\(Q}^\x} + (x\(Q}^\x} + ... (53) 

coincides with the former derivation. The fact that leading finite nuclear mass corrections 
to the matrix elements can be obtained from the adiabatic nuclear equation simplifies their 
derivation. 

For the corrections to the second order matrix element we will need Sq\x) which is 
5M = [E a -H Q -S a (R)-S cl (R)}' [ (0CllWcl) 

Consider now the second order matrix element (Qi Q%) in Eq. ( 1441) with two electronic 
operators Q\ and In order to find the Born-Oppenheimer form and the finite nuclear 
mass corrections, we take the second order derivative 8q 1 q 2 of Eq. (l2Tj) . and multiply from 
the left by (x\ 

(xl<W 2 £a(#) + S QlQ2 S el (R) - 5q 1 q 2 E Su \ x ) + (x\S Ql £ a (R) + S Ql S e i{R) \5 Q2 x) 

+ (X\5 Q2 £,(R) + S Q2 £ el (R) \5 QlX ) = . (55) 

This second order matrix element (Qi Q 2 ) is identified with Sq^E^ and thus 

(Q1Q2) = (4>\Qi {E l _ H)l Qal0> + MQa , E \ H)/ Qi\4>) 

= (x\S Ql Q 2 £c\{R) \x) + (x\ S Qi £ d.( R ) \ 6 Q2X) + (x\Sq 2 £ci{R) \6 Ql x) 

+(x\5 QlQ2 £ a (R)\x) + {x\8qMR)\8qiX) + {x\8qMR)\Sq,x) +■■■ 

= (giQ 2 ) (0) + (Qig 2 ) (1) + ..., (56) 

where Sq-^q^i and 5g 1 Q 2 £ a are corrections to corresponding energies due to electronic op- 
erators Qi and Q 2 , and 

(QiQ 2 > (0) = (x\(MQi (P 1 rr V QM\x) + 

+<xi «w2> [£ ,_ gi ,_ g ' (fi) _ £dWr mf w + - , (57) 
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and 



(QiQ 2 > (1) = <x\(MQi 



In the case (Qi)ei 
subsection. 



(£ e i — H e 



(£ e i — H e \ i 



wQ 2 |^el)|x) 



1 



(£ e l — iJ, 



el 



wQ 2 |0el)|x) 



(£ e i — H e {) 

+ (X\ MQl {F 1 „ v (Q 2 - (Q 2 )^) 1 ff V ^ |0el> |X> 

(£el ~ -"elj (^el — -"el J 

77^ F7 5 77^7 W2 } IX> 



[Sa-fTn-fa^)-^)] 



+ <X| «?1>S 



(0) 



7(Q 2 )^|X)+C.C. 



(58) 



above formulae coincide with those from the previous 



D. Matrix elements of Hermitian operators with /^-derivatives 

Consider an operator of the form Q = Q ■ P where P = — % Vr, and Q is a Hermitian, 
electronic operator, such that [P l , Q l ] = 0, for example the electron-nucleus Breit interaction 
Eq. (1651) . Its matrix element between different rovibrational states is of the form 



(0elXf|Ql0elXi) = g (Xf|[(V/?0el|Q|0el) ~ (0cl|Q| V/J0 c l) J 



+ 2 



(v»KQ).ilxi>-(»KQ).i|Vxi> 

If Q is an imaginary operator, then the second term vanishes and 



(0elXf|Q|0elXi) = (Xf\(Q)e\\Xi), 

(Q)cl = ~i (0cl|<5|V fi e l) = i (V R (f) e l\Q\(j) e l) 

If Q is a real operator, then the first term in Eq. (1591) vanishes and 



(59) 



(60) 
(61) 



kxAQlfaxi) = (xx\(QU-P\xi) = (xt\P-(QUxi)- 



(62) 



IV. RELATIVISTIC RECOIL CORRECTION TO ROVIBRATIONAL ENERGIES 

This is the first and the simplest application of the nonadiabatic perturbation theory, the 
finite nuclear mass correction to the relativistic energy. The total relativistic correction to 
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the binding energy of a E state, neglecting interactions with nuclear spins and the higher 



order 0(m/M) 2 terms, is given by [8 



i, a 



SH = 5H cl + 5H n , (63) 

a,X a>b a>b 



8H e \ = a 2 



(64) 

iH - = *$^4£ + *g)*- («) 

The relativistic correction E Te \ to the energy, taking into account the transformation in Eq. 
(HOD, is 

E vel = (<f>\5H\<f>) (66) 
= (x\(SH)f\x) + (x\(6H)S>\x) + 0(m/Mf\ (67) 



where 



(8H)f = {4> A \5H A \4> A ) = (8H cl ) ch (68) 
{SH)™ = (MSHM-ml^--^-) (0ei|[E^-Vii,^ci]|0ci) 

\1TLA m B / a 

+2 (0 el |5# el 1 ZJ . / # n |0 el >. (69) 



As it have already been noticed in [5j, relativistic recoil effects are of order 0(m/M) and 
can be expressed as a correction (5H)^ to the BO energy E e \(R). 

Since the wave function does not depend on Rg, the momentum Pq implicitly present 
in 5H n can be set to 0. Moreover, we chose Rg in dependence on the particular operator in 
5H C \, in such a way that the result of the transformation cf) in Eq. (ITUl) . namely the second 
term in Eq. ( l69i) vanishes. For example for the first, third and fourth term in 8H e \ Rq is 
the nuclear mass center (ea = + m B ),EB = rn B /(mA + mg)), for 8 3 (r aX ) Rg is 

placed at the nucleus X (e x = 1), and for SH n in the geometrical center [ea = £b = 1/2). 
In this particular choice of Rg, the derivative d5H c \/dR that comes from if^ in Eq. fl69|) 
also vanishes and the relativistic finite nuclear mass correction can be rewritten to the form 

(6H)$ = {fa\6H n \fa) + 2 (MSH el 1 fa + — *0 I0 el ) 

(c G i — ileij V Zm n K / 

1 dV 1 1 dV 

+ ^ {M dR (s cl - HcX y {6Hcl - {6HclU (S el -H el )>* dR M 

1 u \xtt 1 d(V-S cl ) 1 0V |JL 
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The expectation value of SH n is calculated according to Eq. (1321) . namely if 

8H n = Q 1 -V R + Q 2 , (71) 

with Qi electronic operators and VrQi = 0, then 

1 dV i 

(4>el\5H n \(j) el ) = (4> c l\Q 2 \<j) e l) + (<f>d\ft ■ Qljz 7T^7^l^l) + ~= (4>el\n X Ql ■ J n \(j)el)- (72) 

[tel — li e \) Oti K 

One may expect significant cancellations in the leading finite nuclear mass correction (SH)^ 
between the first and the second term in Eq. fl69l) . For example for separate hydrogen atoms 
in the ground state, this correction vanishes. Therefore, the next order correction which is 
0(m/M) 3//2 , may become relatively significant. It is due to, for example, the second order 
in H Q nonadiabatic correction given by Eq. (1371) . or by the orbit-orbit interaction between 
nuclei. 

V. ELECTRIC PROPERTIES 

We will study here the nonadiabatic corrections to the transition dipole moment and 
the electric dipole static polarizability. The direct nonadiabatic calculations have only been 
performed for simple molecules like in jiol . [ill . H . H 2 in 1|, and LiH in [sj. There 
is a lot of literature on electric properties of molecules, in what we call here, the adiabatic 

n 

approximation. Let us mention the extensive review of Bishop [141 and earlier works of 

n n 

Briefer in PJ, 06[. We recover here their results for the electric dipole polarizability, and 
present closed formulae for the nonadiabatic corrections. 

The interaction of a neutral molecular system with the homogenous electric field E is 
given by 

5H = -D- E, (73) 
where the electric dipole operator D is 

D = eJ2x a + e A x A + e B x B , (74) 

a 

and where xa = £b R and x B = —£aR- We note that for charged molecular systems, 
the electric dipole moment has to be defined with respect to the mass center of the total 
molecule. 
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Let us fix the reference frame to the center of the nuclear charge €a = £a/(za + cb), 
e_B = eB/(eA+es), then D = e J2a %a and this D is not affected by the unitary transformation 
in Eq. fllOp . We will show below that interaction of the molecule with the homogenous 
electric field, in spite of including leading finite nuclear mass effects, can be effectively 
described by the Hamiltonian H e Q in the nuclear space 

H cS = -D cl (R) ■ E - E l E> , (75) 

where D e \(R) is given in Eq. (ITHjl and aJ^R) in Eq. (!82|) . From the above equation, the 
transition dipole moment is 

D & =<Xf\D*\Xi), (76) 

and one notes that the matrix elements between the same nuclear states always vanish. 
Similarly, the total electric dipole static polarizability is 

^ = (xKilx) - 2 (x\D l cl 1 - - DM, (77) 

where the first term is due to the electron excitations and the second one is due to the 
rovibrational excitations. For molecules, with so called permanent electric dipole moment, 
the second term is dominating, particularly significant contribution comes from intermediate 
nuclear states with the same vibrational number, but different J. 

The electronic matrix elements D e \ and a e \ are obtained as follows. According to Eqs. 
( [351401 and Eqs. ( 12911301 the electric dipole moment including the leading nonadiabatic 
effects is 

D fA = D$ ) +D%\ (78) 

where 

= (0cl|^|0el), (79) 

= 2 (MD (f 1 „ v H n |0 el ). (80) 
For E electronic states the nonadiabatic correction can be rewritten to the form 

(£ei - H cl )' \ 2m n R I J m n R 1 (£ el - H A )' 

j-u in 1 d ( v ~^) 1 9V 
m n {(pclllJ (£ el - H el )' dR (£ cl - H c ^ dR l</W 
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This result is in agreement with the previously obtained much simpler formula for the HD 
molecule where only the last term in H n contributes in the above equation, due to the 
inversion symmetry of the ground electronic state. We observe that in spite of including first 
order nonadiabatic corrections, the transition dipole moment can be represented in terms 
of the matrix element of the electronic dipole moment _D e i evaluated with nuclear functions 
Xf and x\i see Eq. (|39|) . A similar result holds for the electric dipole static polarizability 



ij (0)ij , 

«el = "el +"cl , 



where in the BO approximations it is 

-2(0 e l|^ 



(0) ij 



(Sol - #el)' 

and the nonadiabatic correction, using Eq. (|58|) is 



Wei), 



(82) 



(83) 
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«el 



5el 



D l 
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(£ c i — H c \ 

+ (0el|#n ,„ 1 TT , l (D i -D™ i 



^|0el) 



el 



+ ( 



D l 



(£ e i — H d 
1 



-{D 3 — D x 



'el 



(0)^ 
el i 



(£el — 
1 



J7^l0el) 
H n |0 e l) 



(84) 



(£ e l — H e i)' 6 (£ c i — H e \) 

The explicit formula after taking the derivative with respect to R is little bit too long to 
be written here, but can easily be obtained with the help of Eqs. ( l2"Uf30ll3ip . Eq. f|77|) in 
the adiabatic approximation, namely a e \ — > , D e \ —>■ was already obtained in the 
literature Results for nonadiabatic corrections when a e \ — > oi^\ D c \ — > are, 

in our opinion new and have not been discussed in the literature. 



VI. MAGNETIC PROPERTIES 



We consider here nonadiabatic corrections to magnetic properties of a diatomic molecule. 

n 

These properties within the BO approximation are reviewed in details by Flygare in [18]]. 
Here, we demonstrate that the nonadiabatic corrections can be implemented in the effective 
nuclear Hamiltonian, similarly to that in the BO approximation, namely 

H cS = - V (R) I- J — 7/ 1(1 - a el (R))B - lJ {R) J B - ~ B X (R)B, (85) 

where it is understood that the electronic part of the angular momentum operator J e \ Eq. 
ffT7|) vanishes on x, V is the spin-rotation constant, fij = 1 7/ is the nuclear magnetic moment, 
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[ij = J is the orbital magnetic moment, a c \ is the independent shielding constant, \ is the 
magnetic susceptibility, and B is the magnetic field. The Hamiltonian in Eq. fl85l) should in 
principle involve also quadrupolar interaction (P Pp 2 ' (J* Ji)^ or (P Pp 2 ' J % BK The first 
term comes from the quadrupole moment of nucleus or tensor interaction between nuclear 
magnetic moments, while the second term has not been investigated in the literature so far. 
Both terms will not be considered in this work. If two different nuclei are involved, each one 
has its own spin and magnetic moments, then interaction between them should be included, 
but we do not consider it either. In the following we rederive in a simple way known results 
for rotational magnetic moment /ij and spin rotation constant 77 and obtain nonadiabatic 
corrections to the shielding constant. These corrections, up to our knowledge, have not yet 
been investigated in the literature 



A. Nonrelativistic Hamiltonian in the external magnetic field 

In order to obtain finite nuclear mass corrections we start with the Hamiltonian of a 
molecular system in the homogenous magnetic field 



#o = £^-+V, (86) 



where f3 enumerates electrons and nuclei. We assume the coordinate system as defined in 
Eqs. (EJHj): X/3 = fp — Rq and perform the following unitary transformation 

H = e-^H e^ + d t( p, (87) 
<f = J2 e P duxp- A(R G + uxp) 

(3 J ° 

= A* + 44^/2], (88) 

13 

= + (89) 

where A = A(Rg) ■ The result of this transformation on the canonical momentum is 

e-^Ke^ = P i + ^(S a xBy, (90) 
e-^vr^ = p ^ + ^(x A xBy + e -^(DxBy, (91) 
e-^e** = Pb + y @b x By + 6 f (D x By, (92) 
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where D is the total electric dipole moment D = e /3 %p an d we assumed that J2a e p — 0) 
so the molecule is neutral. In consequence, the transformed Hamiltonian does not depend 
on Rq, thus we are allowed to set Pq = 0. In order to further simplify, we perform the 
next transformation as defined in Eq. (flOl) . neglect 0(m/M) 2 terms, and the transformed 
Hamiltonian takes the form 

H = H cl + H n + H IM + H x , (93) 
Hp = -J2 T^^aXQa-B + SH^ (94) 

SH, = -(I±-1^)dxP-B-U— + —) (e A e 2 B + e B e 2 A )RxP.B 
\m A m B J 2 \m A m B J 



+ I(A + A)fl xEa .a + 5^(£i_f!L)a xI ; 4 .a i ( 9 5) 

2 \m A m B J a 2 \ m A m B J ^ 



c 2 



H x = E (*« x B) 2 + 6H X , (96) 



a ^ m 



SH X = — (e A x A x B + e A D x B) 2 + — !— (e B x B x B + e B D x B) 2 . (97) 
8mA 8m B 

The center of the reference frame Rq is placed arbitrarily on the symmetry axis. This 
freedom will be used below to simplify formulae for nonadiabatic corrections. 



B. Rotational magnetic moment 



The rotational magnetic moment has been first investigated by Wick in [20| for H2, later 
by Ramsey in who extended Wick's result to arbitrary masses of nuclei and performed 
first measurements in Ref. 2l| and improved measurements with Harrick in Ref. 22] . 
The rotational magnetic moments results from coupling of the molecular rotation to the 
magnetic field. The expectation value on the E state of the first term in Eq. (j94p e/(2 m) J e \ ■ 
B vanishes, thus the leading coupling comes from nonadiabatic corrections to the matrix 
element as given by Eq. (j4"3l and from SH^, namely 

7j = (— ~ —) + ^ ^ + &B — - — 1 (0 el | J e i 1 j d U el \ ( 98 ) 

\m A m B J R 2 2m n 2mm n R 2 e e (£ e \ — H e i)' e 6 

In the center of nuclear mass frame e A = m A / (m A +m B ), e B = m B /(m A +m B ) the rotational 

magnetic moment in units of the nuclear magneton /i/ becomes 

1j fJ-j fZ A Z B \ m p 

= = rn m -4- — 



„, ~ J H ~ m " m ' \^A + mi) + l^W {MJ '< Je^H^ J « M < (99) 
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where m p is the proton mass, in agreement with the result from Ref. [21] [Eq. (4) with 
assuming Za = Zb = 1]. One observes that rotational magnetic moment in Eq. fl98|) does 
not depend on the choice of reference frame, but depends on the distance R between nuclei. 
For large R it vanishes at least with R~ 6 . The dependence of jj on R leads to the appearance 
of the magnetic dipole transition in if 2 , between states of the same angular momentum, but 
of different vibrational number. These transitions, up to our knowledge, have not yet been 
investigated in the literature and may play a role in the astrophysical environment, however 
their importance should be verified by explicit calculations. 

As was first noted in Ref. 22J the rotational magnetic moment is related to the param- 
agnetic part of magnetic susceptibility x- m the BO approximation \e\ is given by 

xS = ~ E <M#> - «l« - (MJi (ioo) 

For E states it can be simplified to the form 

Xei = ^(^l-in-Xafl^n'^ + i^xl + in 

8m „ 

» 2 i 



-I^^ J >'Ie^Hj> J >^'> l - s ''- nW) ' (101) 

The last terms in both Eqs. (|98|) and (11011) . which are the second order matrix elements, 
are similar while the first terms in both equations are simple to evaluate. This allows one to 
express the difficult to measure magnetic susceptibility in terms of the rotational magnetic 
moment. This relation however, works only in the BO approximation, as nonadiabatic 
corrections will be different. As noticed by Authors of [23|], this second order matrix element 
with J e \ operator, is present also in the nonadiabatic equation for the nuclear function x as 
the W± function (in the notation from our previous work 



C. Spin-rotation Hamiltonian 



The general spin-orbit Hamiltonian for arbitrary nuclei, including the external magnetic 
field is (h = c = 1) 0, S Q 

e a ep 1 



5h=y: 



s a ■ r a/ 3 x -rip — - — s a ■ r aP x n a 



(102) 
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where summation over a, (3 goes over electrons and nuclei. In particular, the coupling of 
the nuclear spin I — sa to the rotation and to the magnetic field using Eq. (I102p is 

(9a - 



5H 



E 



e A e 



b 4vr 2r\ b 
e A e B I 



4tt 2r AB 



9A 

itla m 

9a 
m A m B 



TAb X 7T b 



tab x ttb 



m , 



(9, 



TAb X 71 A 
-I) . 



m 



tab x 7i A 



A 



(103) 



For convenience we chose the reference frame centered at the considered nucleus Rg = Ra, 
so €a = 1,c_b = 0, perform unitary transformations in Eqs. (I88|ll0p . and obtain 



5H 



E 

b 

eA e B 



e A e 



47r 2 m a 



%b 



Qa ( ^ m - e 
— (% + P + -x b xB 



R 



x 



4vr 2m A R 3 



m 

9a 
m B 



P + — RxB 
2 



(9A - 1; 



(gA - 1) 



771,4 



P - g e i 

D 



x B 



(104) 



where = J2 a Qa, %e\ — J2a%a, and the electric dipole operator is D = ex e \ — e B R. This 
Hamiltonian will be used in next sections to rederive the known formulae for the spin- 
rotation and the shielding constants, and to obtain a closed expression for the nonadiabatic 
corrections to the magnetic shielding, which contribute at the level of m e /m n what for H2 
is about 10 -3 . 



D. Spin-rotation constant 



The theory of the spin- rotation interaction was introduced by Wick 20j, Ramsey 21], 
and Foley 25j, and further developed by Ramsey 2y, |27|, by Frosch and Foley [28J and 
again by Ramsey 2£|. These theoretical results were not in agreement for a long time, until 
Reid and Chu in 30j found a complete set of corrections. Here we rederive their result on 
the basis of Eq. (TT04|) . 

The spin-rotation interactions results from 5H above with B = 0. For considered S states 
terms with q e \ vanish and 6H takes the form 



5H 

Qi 
Q2 



Qi ■ I + Q 2 x P ■ I, 

eA_e g A x b x q b 
^ 4 7r 2mm a x 



eA e 1 Xb 



, 4tt 2m\x\ 



E 



6 



1 



4 7T 2)71,4 



_9^4_ (gA - 

m B m A 



R 
W 



(105) 
(106) 

(107) 
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The expectation value of (0 e i|Qi|0ei) vanishes, and Qi operator contributes only through 
nonadiabatic matrix elements Eq. (1431) . so the total spin- rotation constant is given by 

J- J , . . * 1 



r]I-J 



m n R 2 



^el|<^el 



r.QM + MQM xP-i. (108) 



(£ oi — H c \ ) 

The expectation value of the first term in Q 2 can be expressed in terms of derivative of BO 
energy, namely 

fdS e i e A e B 1 



1 3 ^1 

b 47F X b 



11 



\dR 4vr R 2 
and thus rj in atomic units becomes [ex = —Zx e, a = e 2 / (47r)] 



(109) 



1 Z A 9 A 

R 2 2m n mA 



J2 X 9a 



1 ^ f a X g b | 

(£ C 1 — -f^el)' 5 



X) 



l) 



1 1 



d£ e \ 1 Z A Z B g A 



R2w 2 A dR ' 2m im J' 
in agreement with Ref. 30[], their C is related to our rj by C = 2 7rm Since, there is 
cancellation between the first and the third term, r\ vanishes at least as R~ 6 for large values 
of R. 



E. Magnetic shielding constant 



The shielding of the external magnetic field due to atomic electrons, was first considered 



by Ramsey in 



26j with the help of the Breit-Pauli Hamiltonian with the external magnetic 



field. Here we rederive his result and in addition to it, obtain nonadiabatic corrections. A 
sim ilar calculations for ato mS have recently bee, performed in @,Q. 

The Hamiltonian of a molecule in the magnetic field including the nuclear spin, but 
neglecting that of electrons, is a sum of H in Eq. (193]) and SH in Eq. (11031) . In the BO 
approximation it? is a sum of diamagnetic and paramagnetic parts. In atomic units, they 



are correspondingly 



T»fB 



a 



9el 



X b 



Xi 



{x b x B)\(j) el ) 



+ (4> c i\J2x a x q a - B 



{£ — H t 



el J 



Xb x q b 



x b 



111) 



Nonadiabatic corrections , namely all corrections which are linear in the electron-nuclear 
mass ratio come from several sources, and we split them into four parts 



(i) 

cl 



r (l) 

n 



(i) 



;ii2) 
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o~W is the correction due to H n to the matrix elements in Eq. fillip , namely using Eq. 
one obtains 



Iai 1] B = a 2 



(S - H el )' 



n <Pc\ 



(£-# e 



■I'l, 



+ (0el| #n 



./, iv-, - 5 1 v^ X * r 1 



(£-# e 



" #el)' 



n Pel/ 



(113) 



a - ,} is a correction to the diamagnetic part due to the direct coupling of the nuclear spin to 
the magnetic field in 5H, Eq. (11041) 
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<> el >. (114) 



<3\P is a correction to paramagnetic part due to spin-rotation interaction in 5H Eq. fjl04j) 



6 X 6 



P + (<M - 1) q c 



m B m A 



(115) 



The derivative over nuclear coordinates P in the above and below, can act on the right or 
on the left, since these matrix elements does not depend on this. Finally, Oj is a correction 
to the paramagnetic part due to SH^ in Eq. fl95l) . 



Ta[ l) B = a 1 



Ml— (xci + Z B R) x (2P - q cl ) - ^Rx P \ ■ B 



■T, 



(116) 



The total magnetic shielding a is obtained by averaging with the nuclear wave function 



117) 



and one notes, that the orbital magnetic moment with the spin-rotation coupling gives 
additional contribution to the shielding constant 



Sa» = ~(X\V J 1 „ 1 F FV 1l J j \x), 

7/ [Eg, - H n - £ a - 8 C \\' 



(118) 
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which however is negligible. 

Let us present in more details the averaged shielding a = cx n /3 in the case of H2 and 
isotopomers Z A = Z B = 1, 

*=(X\<7$ ) +4 ) \X), 

and 
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(120) 



(121) 
(122) 



(123) 



(£ -H t 



7 2^ -3 lw- 



(124) 



These formulae can be further simplified by shifting the reference frame to the geometrical 
center and by using gerade symmetry of the ground electronic state of H 2 and isotopomers. 



VII. SUMMARY 



We have presented a general approach to finite nuclear mass corrections in molecular 
properties which is based on the nonadiabatic perturbation theory [3j. These corrections 
were represented in terms of electronic matrix elements averaged with the nuclear wave 
function, similarly to that in the adiabatic approximation. We obtained formulae for nona- 
diabatic relativistic corrections which can be used to perform accurate calculations of dis- 
sociation and rovibrational energies. Currently the accuracy of theoretical predictions in 
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H2 \^ are limited by these not well known effects. Similarly, we obtained formulae for the 
nonadiabatic corrections to the transition electric dipole moment and the electric dipole 
polarizability. They can be used for the comparison with precise measurements of polariz- 
abilities, for example in such complicated system as excited vibrational states of the water 
molecule I3J 



]. Finally, we presented nonadiabatic corrections to the magnetic shielding, 
which are important for molecules involving hydrogen or deuterium, where the finite nuclear 
mass significantly ~ 10~ 3 affects the magnetic shielding. 
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